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Abstract. We report on a comprehensive systematics of fusion-evaporation and/or fusion- 
fission cross sections for a very large variety of systems over an energy range 4A-155A MeV. 
Scaled by the reaction cross sections, fusion cross sections do not show a universal behavior 
valid for all systems although a high degree of correlation is present when data are ordered 
by the system mass asymmetry. For the rather light and close to mass-symmetric systems the 
main characteristics of the complete and incomplete fusion excitation functions can be precisely 
determined. Despite an evident lack of data above 15A MeV for all heavy systems the available 
data suggests that geometrical effects could explain the persistence of incomplete fusion at 
incident energies as high as 155A MeV. 



1. Introduction 

Three regimes of fusion process only qualitatively explain the fusion excitation function [1]. 
The two low energy regimes have been thoroughly studied in the past as well as today [2]. On 
the contrary, the third high-energy regime related to the extinction of fusion cross section has 
been studied only occasionally. The question how and why the fusion process disappears did 
not receive a satisfactory answer yet. The two systematics on fusion excitation function data 
have been published so far, in 1984 by Mogenstern et al. [3] and by the Indra collaboration in 
2006 [4] both being focused on rather light and symmetric systems. In the former study, the 
disappearance of complete fusion has been foreseen for vl < 0.19c and the onset of incomplete 
fusion for 0.06c <V[}- For a symmetric system and in terms of incident energy these limits are 
respectively equal to Q7A MeV and 6,7^4 MeV. In the present work we challenge both previously 
established limiting values. In the latter systematics, presented is a fusion excitation function 
for 7 light nearly symmetric systems. It has been found that quasi-fusion becomes vanishingly 
small above 50^4 MeV. 

In this work we have considerably extended the previously published systematics and an 
exhaustive review of measured fusion cross sections available today at energies well above the 
Coulomb barrier (-Slab > 4^4 MeV) is presented. Undertaken a scrutiny scan of the published 
fusion data over the last forty years ended with 168 fusion cross section (FCS) values belonging 
to 57 different systems: A span in total system mass ^4tot = At + A p {At and A p stand for 
target and projectile mass, respectively) ranges between 26 and 246, system mass asymmetry 



1 vl denotes the velocity of the lighter partner in the center-of-mass frame. 



ji = \ At — A p \/(At + A p ) from to 0.886, neutron to proton N/Z ratio from 1 to 1.536 while 
incident energy lies between AA and 155^4 MeV. 

Created high-energy excitation function for the fusion-evaporation and/or fusion-fission 
processes should allow 1) to pin down a possible regularity in the global available data set 
or in appropriately selected data subsets, 2) to identify a possible missing data and to establish 
a priority list of new measurements that would have to be performed and 3) to deliver a set of 
data that can be used to constrain microscopic transport models used to describe properties of 
heavy-ion collisions around the Fermi energy. In fact, reproducing the overall evolution as well 
as the main features (onset, maximum and extinction) of the global fusion excitation function 
or the one obtained by sorting data as a function of a given entrance channel parameter is a real 
challenge and a sine qua non condition for these models. Indeed, before pretending to describe 
in detail properties of the reaction exit channel such as creation and properties of particle species 
with their energy spectra and isospin these models firstly have to accurately reflect the global 
features of reaction mechanisms involved. 



2. The raw data 

In Fig. 1 are displayed all collected fusion data as a function of incident energy in the laboratory 
reference frame sorted by the increasing system mass Atot- The systems are differentiated 
by a color code (online): blue and green symbols label the lighter systems (26 < A tot < 116) 
while pink and red symbols the heavier ones (146 < Atot < 246). Black symbols are used to 
point out measurements for which in the original publication is explicitly mentioned a possible 
overestimation of the fission component due to an unresolved contribution of fast-fission [5— 
8]. Black symbols are also used for the five systems labeled by Only ER in Fig. 1 for which 
the evaporation cross sections were measured only, whereas a more or less significant fission 
component was expected too [9, 10]. 
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Figure 1. (Color online.) Raw 
fusion cross sections plotted as 
a function of the laboratory en- 
ergy per nucleon E^\b/A. The 
inventoried systems are distin- 
guished among them by a color 
code which is used throughout 
this paper. For convenience, all 
the references are indicated in 
Figs. 3 and 4. See text for de- 
tails. 



It is instructive to investigate the systems by sorting them following both their mass 
asymmetry and their isospin. By a careful survey of the list of systems one infers easily that 
lighter systems (top of the list) are rather symmetric both in mass and isospin (0<fi<0.5 and 
1 <N/Z < 1.25) while, on the contrary, heavier systems (bottom of the list) are very asymmetric 
in mass (for most of them /i>0.75) and, as expected, more asymmetric in isospin (N/Z> 1.3). 
The ensemble of fusion data splits into the two distinct sets exactly matching our color code. A 
first set corresponds to Light Symmetric Systems (LSS) for which FCS diminish with incident 
energy and disappear above 40-50^4 MeV (see also [4]). A second set corresponds to Heavy 
Asymmetric Systems (HAS) for which FCS increase up to 20A MeV and then decrease. The 
HAS branch surprisingly persists at energies as high as 155A MeV. At that point one must 
stress an evident lack of FCS data: For the heaviest systems in general data are scarce above 
15A MeV, particularly for HAS between 30^4 and WOA MeV, and for systems of medium mass 
asymmetries in the entire energy range. Let us add that fusion data above 100A MeV are from 
a single experiment [8]. 
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Figure 2. (Color online.) 
Kox over Tripathi reaction cross 

The 



Kox /^Tripathi 
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74 



section a 
black vertical line indicates the 
Ec.m./A =1.6A MeV limit. See 
text for details. 



3. Normalization of the fusion cross sections 

To reduce the effect of Atot one normalizes FCS by reaction cross sections o"f us /crR. By specifying 
abscissae in units of available energy (i.e. the system center-of-mass energy per nucleon E c . m JA) 
one expresses the mass asymmetric systems on the same footing with those which are mass 
symmetric. These coordinate scalings greatly ease the comparison of various systems with each 
other. Several parameterizations are used in calculation of reaction cross sections [11-16]. The 
parameterization of Tripathi is chosen because it is applicable for any combination of colliding 
nuclei and valid over the large energy range from a few A MeV to a few A GeV [11]. Note that 



this parameterization depends on the root-mean-square charge radii of nuclei for which we used 
the tabulated values of Ref. [17]. Nevertheless, to check that our results do not depend too much 
on the chosen prescription, a comparison between the Kox [12] and Tripathi formulas has been 
performed, Kox formula being modified in order to extend its applicability down to a few A MeV 
similarly as Shen et al. have done [13]. The comparison is displayed in Fig. 2 where is shown 
the ratio a^ ox / a~j^ npathl of reaction cross sections calculated respectively by Kox and Tripathi 
formulas. Again, the two sets previously mentioned (LSS and HAS) are evident. For HAS, it is 
clear that this ratio is almost constant, varying between 1.09 and 1.17 for E c . m JA> 1.6A MeV. 
Below this limit denoted by the thick vertical line in Fig. 2, the differences are growing and 
reach more than 60 % at the lowest energies. 

For LSS, the conclusion is similar except that the ratio is closer to 1 (within 10%). 
Subsequently, we will restrict the domain of validity of our analysis to E c _ m _/A > 1.6A MeV 
where both formulas are quite compatible. Thus, choosing Tripathi formula as we have done, 
implies a slight difference in absolute values of normalized cross sections but does not alter the 
shape of the excitation function. Anyway, it is likely that none of these parameterizations is 
correct at such low energies. To recall this zone of inconsistency of evaluated total reaction cross 
sections <tr in all next figures a yellow band is added for E c _ m JA< 1.6A MeV. 

Figures 3 and 4 show the LSS and HAS data sets in the scaled coordinates, respectively. A 
particularly huge effect is observed for HAS, but clearly, contrary to what one may expect, FCS 
cannot be put under a universal behavior that would be valid for all systems. In fact, instead 
of following a single global law, plotted FCS strongly suggests the existence of two branches: 
a lighter symmetric one (Fig. 3) and a heavier asymmetric one (Fig. 4). These figures call for 
several comments. When only LSS are considered (Fig. 3), apart from a few points corresponding 
to the N + C and S + C systems, a very regular correlation is observed. The preliminary FCS 
data on a symmetric and rather heavy Xe + Sn system reported at this Conference [62] seems to 
nicely fall in the mentioned systematics (not included in Fig. 3). The LSS excitation function 
can be very well reproduced by a hyperbolic function (see the solid red curve in Fig. 3). 

The fit was achieved with the points located outside the yellow zone. The fit parameters, 
however, are quasi-independent on whether or not the points of the yellow zone are taken into 
account. The result is also independent on whether the Kox or Tripathi formula was used. Both 
<7f us and ctr depend on available energy. No obvious and simple physical explanation justifying 
the observed hyperbolic dependence can be drawn using e.g. 1/E c . m , dependence as suggested in 
the Bass model [63] . It is probably due to a non- validity of this model at the highest energies of 
the so-called regime III of fusion process (see e.g. in Ref. [1]). Further important result is that 
incomplete fusion reaction mechanism clearly disappears around i£ c .m. = 12-A MeV whatever the 
LSS system is. This findings is in agreement with the results of Ref. [4]. By using the microscopic 
transport Landau- Vlasov model [64] we have shown years ago that the fusion process abruptly 
disappears around the Fermi energy in head-on collisions: An elongated composite system gives 
rise to two nuclei of the reaction exit channel having kept a strong memory of the entrance 
channel suggesting a partial transparency phenomenon [65]. 

Understanding the HAS data subset seems to be less straightforward. Nevertheless, by 
neglecting the data points of the yellow zone and those labeled by Only ER the few remaining 
data suggest the existence of a second asymmetric branch of FCS tending towards a constant 
value at high incident energies. Can one explain the persistence of fusion at such high energies? 

Above 100^4 MeV a crucial role that the reaction geometry plays in heavy ion collisions 
is a well established fact. In such geometrical picture, fusion can occur for collisions with 
an impact parameter smaller than 6 maX i a value corresponding to a complete overlap between 
the projectile and the target. Otherwise, for impact parameters greater than 6 max , part of the 
projectile pursues its trajectory at small angle with a velocity close to the projectile one. Within 
this approximation and with a reaction in direct kinematics the fusion cross section is expressed 
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Figure 3. (Color online.) Normalized fusion 
cross sections for lighter and rather symmetric 
systems plotted as a function of E c , m ,/A. See 
text for details. 



Figure 4. (Color online.) Same as Fig. 3 
for heavier asymmetric systems. See text for 
details. 



as following 



vr hi 



7T {Rt - R t 
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(1) 



where R p (Rt) is projectile (target) radius. Assuming the simple geometrical formula for the 
reaction cross section 

a R = vrr 2 + A 1 / 3 ), (2) 



one gets 
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■ 1— a I 



(3) 



Eq. (3) indicates that at high energies normalized FCS tend towards a value which only depends 
on the mass asymmetry \x. For the 14 N + 154 Sm and 14 N + 197 Au Eq. (3) gives 14.4 % and 
17.0 %, respectively in the very good agreement with experimental data. This reaction picture 
is confirmed by our Landau- Vlasov simulations [66] : For impact parameters of complete overlap 
collisions lead to the formation of a massive incomplete fusion nucleus accompanied by pre- 
equilibrium emission of light particles in the forward direction, an emission mostly originating 
from the projectile. Obviously, for a more and more symmetric system this limit tends to zero in 
accordance with what is observed experimentally. The mass asymmetry fi of the system seems 
to be the key parameter in explaining the above transition, at least in the high energy part. The 
data on systems of medium \x are needed to confirm or deny the above assumption. 



4. Complete fusion excitation function 

Among the experiments listed in our systematics only ten have explicitly been designed to 
measure both complete and incomplete fusion components [21, 25-27, 30, 31, 35, 36, 38, 43]. This 
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Figure 5. (Color online.) 
Normalized complete fusion cross 
sections for lighter symmetric 
systems plotted function of 
the available energy -Ecm. /A. See 
text for details. 



allows establishing a complete fusion excitation function incorporating 29 FCS data points 
relative to 10 systems spanning a total mass range from 40 to 68. Data are reported in Fig. 
5. Again, a very regular correlation is observed with FCS decreasing rapidly as a function of 
Ec.mJA. As before, it can be fitted by the same kind of hyperbolic law (see pink solid curve 
in Fig. 5) which crosses abscissa around QA MeV. For the complete fusion reaction mechanism, 
center-of-mass energy and excitation energy are equivalent. Thus, the energy limit of 6 A MeV 
corresponds to the maximal excitation energy which can be deposited into a light compound 
nucleus. The above conclusion rises several questions. What happens for heavier systems? What 
for asymmetric systems? Does one obtain the same behavior and the same limit for complete 
fusion independently of Atot? Again, further measurements would be welcome. 

Finally, in Fig. 5 are superimposed the two hyperbolic fits obtained for complete fusion (CF) 
and for the sum of CF and incomplete fusion (IF) components and which are denoted by pink 
and red curves, respectively as well as their difference which is colored blue. These functions 
provide a measure of the average contribution of CF and IF fusion components and of their 
evolution with energy. From the results of Fig. 5 one can, in particular, infer 

• The threshold of IF is around E c , m JA equal to 1A MeV. 

• By increasing E c _ m JA CF decrease while IF increase. At AA MeV each of these two 
components cjif and ocf contribute about the same value of ~ 0.1 x cr to the total 
reaction cross section. 

• IF component increases up to QA MeV where it reaches a maximum while at the same time 
the CF process vanishes, ajp reaches 15% of cjr. 

• After this maximum IF decreases and disappears around 12^4 MeV. 



To our knowledge, it is the first time that such accurate limits concerning both CF and IF 



fusion components are deducted from such a large body of data. It allows to draw a simple 
picture for the evolution of these reaction mechanisms. 



5. Conclusions 

In summary, we studied the energy dependence of fusion cross section aiming at understanding 
the evolution of the underlying reaction mechanisms as a function of entrance channel parameters 
and the causes of its (non-)disappearance. The heart of this work is a systematics and as wide as 
possible overview of measured fusion-evaporation and/or fusion- fission cross-sections including 
summed up complete and incomplete fusion contributions. Normalized by reaction cross sections 
and plotted as a function of the available energy per nucleon, the fusion cross sections show a 
rather universal behavior leading to a disappearance of the fusion around 12^4 MeV in the case 
of light symmetric systems. In the case of heavy asymmetric systems incomplete fusion seems 
to persist and tends towards a constant value depending only on the mass asymmetry of the 
system. In addition, dissociating available complete and incomplete fusion cross sections allows 
to provide fairly accurate energy limits characterizing these two mechanisms. Both the observed 
trends of fusion excitation function and the above limits place important constraints on the 
basic ingredients of microscopic transport models. 
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